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Broad-Band X-Ray Study of a Transient Pulsar RX J0059.2-7138 
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Abstract 

We report on the results of the ASCA and ROSAT observations on RX J0059.2— 7138, a transient X-ray 
pulsar in the Small Magellanic Cloud. The barycentric pulse period has been precisely determined to be 
2.763221 ± 0.000004 s. The pulse shape is almost identical in all of the energy bands. The pulse fraction 
increases with the photon energy below ~ 2 keV, while it is nearly constant at ~ 37% above ~ 2 keV. The 
X-ray spectrum has been found to consist of two components. One is dominant above 2 keV, and exhibits 
sinusoidal pulsations. This component is well described by a typical model found in many X-ray binary 
pulsars, a power-law of photon index 0.4 with an exponential cut-off at 6.5 keV. The other is dominant below 
1 keV and shows no significant pulsation. This component is represented by either a broken power-law with 
photon indices of 2.6 and 5.1 below and above a break energy of 0.9 keV, or a metal-poor thin-thermal 
plasma with a temperature of 0.37 keV. The phase-averaged luminosity is ~ 10 38 erg s _1 (0.1-10.0 keV) 
for both components. A hint of oxygen over-abundance is found in the absorbing column, possibly due to 
circumstellar gas ejected from an evolved companion. 

Key words: circumstellar matter — pulsars: individual (RX J0059.2— 7138) — stars: Be — stars: 
neutron — X-rays: binaries 



1. Introduction 

Accretion-powered X-ray pulsars are binary systems 
consisting of a neutron star and a stellar companion 
(here, X-ray binary pulsars or XBPs in short). The grav- 
itational energy of accreting matter is converted to X-ray 
radiation, hence its luminosity is generally variable, de- 
pending on the mass-accretion rate. Transient behavior 
is also observed from many XBPs, most of which have a 
Be star companion with an eccentric orbit (e.g. Stella et 
al. 1982; Bildsten et al. 1997). 

XBPs have been observed mainly in the hard X-ray 
band (~ 2-40 keV) with non-imaging satellites (e.g. Na- 
gase 1989). Their spectra are generally described by a 
power-law with a high-energy exponential cut-off (here- 
after, ECUT power-law) with a photon index of T ~ 1 
and a cut-off energy at around 10-20 keV. A fluorescent 
emission line at 6.4 keV from neutral iron atoms has been 
observed in many bright XBPs, which serves as a probe 
of the circumstellar matter of the binary system. In the 
soft X-ray band (below ~ 2 keV), however, their spec- 
tra have not been well studied due to a large interstellar 
absorption, because most of the XBPs are located in the 
galactic plane (e.g. Bildsten et al. 1997). 

In order to study soft X-rays from XBPs, therefore, the 
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Magellanic Cloud sources have great advantages, owing 
to the low interstellar absorption and well-calibrated dis- 
tances. Woo et al. (1995, 1996) made broad-band spec- 
troscopic studies in the 0.1-40 keV band of SMC X-l and 
LMC X-4, respectively, and found a "soft excess" below 
~ 2 keV, in addition to the usual hard spectrum given by 
an ECUT power-law model. Such a soft excess was also 
discovered from a transient XBP XTE J0111. 2-7317 in 
the Small Magellanic Cloud (SMC), by an ASCA obser- 
vation in the 0.5-10.0 keV band (Yokogawa et al. 2000). 
However, since such a broad-band spectroscopy has been 
performed on only a few XBPs, whether or not the soft 
component is common among XBPs is still unclear, es- 
pecially for transient XBPs. 

Hughes (1994) serendipitously discovered a new tran- 
sient XBP, RX J0059.2-7138, in the SMC with the soft 
X-ray band observation of ROSAT, and found that the 
spectrum is unusually soft, which is composed of a black- 
body with a temperature of kT ~ 35 eV and a steep 
power-law with a photon index of Y ~ 2.4. A proposed 
candidate for the optical counterpart was revealed to be 
a Be star by Southwell and Charles (1996). 

In order to examine the nature of the unusually 
soft spectrum of RX J0059.2— 7138, observations in the 
higher energy band are essential. RX J0059.2— 7138 was 
also detected in a simultaneous observation of the SNR 
E 0102-72.2 with the ASCA satellite, which is sensitive 
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in the ~ 0.7-10.0 keV band. A preliminary short report 
on the ASCA results can be found in Kylafis (1996). 

In this paper, we combine the ROSAT and ASCA 
data, and perform a broad-band timing spectroscopy, 
covering ~ 0.1-10.0 keV. We show that the spectrum 
of RX J0059.2— 7138 has a hard component, which re- 
sembles the spectra of usual XBPs well, and that the 
pulsations are due to the hard component. In addition, 
a non-pulsating soft component is present in the spec- 
trum. We assume the distance to RX J0059.2— 7138 to 
be 60 kpc, the nominal value to the SMC (Mathewson 
1985). 

2. Observations 

The transient XBP RX J0059. 2-7138 was serendipi- 
tously detected on May 12-13 in 1993, in the simulta- 
neous ASCA and ROSAT observations of E 0102-72.2, 
the brightest supernova remnant in the SMC. ASCA 
and ROSAT observations spanned about 99 ks (MJD 
49119. 38-MJD 49120.52) and 20 ks (MJD 49119.95-MJD 
49120.18), respectively. 

ASCA carries four X-ray Telescopes (XRT) with two 
Gas Imaging Spectrometers (CIS 2 and GIS 3) and two 
Solid-state Imaging Spectrometers (SIS and SIS 1) on 
the focal planes (Tanaka et al. 1994; Serlemitsos et al. 
1995; Ohashi et al. 1996; Burke et al. 1994). Since 
RX J0059. 2-7138 was outside of the SISs' field of view 
and was located near the calibration source of GIS 3, 
we used only the GIS 2 data in this paper. The GISs 
were operated in the normal PH mode with a time res- 
olution of 0.0625 s (high bit rate) or 0.5 s (medium bit 
rate). We rejected any data obtained when the satellite 
was in the SAA (South Atlantic Anomaly) region, when 
the elevation angle of the Earth limb was less than 5° or 
when the cut-off rigidity was less than 4 GV, and also re- 
moved particle events using the rise-time discrimination 
method. The total exposure time of GIS 2 was ~35 ks 
after screening. 

ROSAT carries a soft X-ray telescope with a High Res- 
olution Imager (HRI) and a Position Sensitive Propor- 
tional Counter (PSPC) with one of the two on the focal 
plane (Trumpcr 1983). In this observation, PSPC-B was 
on the focus and operated in the pointing mode. We took 
the screened data from the HEASARC archive, of which 
the total exposure time was ~ 5 ks. 

3. Analyses and Results 

For the GIS data, X-ray photons were collected from an 
elliptical region around RX J0059.2-7138 (figure 1). The 
background data were extracted from off-source areas on 
the edge of GIS 2 at the same off-axis angle as the source 
(figure 1), because both the diffuse X-rays and non- X-ray 
background depend on the off-axis angle. 
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Fig. 1.. ASCA GIS 2 image. The sources detected at the center 
and the edge arc the SNR E 0102-72.2 and RX J0059.2-7138, 
respectively. The photons of RX J0059.2— 7138 were spreaded 
by the point-spread function at the detector edge region. The 
source and background regions used in the analyses are enclosed 
by the dashed and solid lines, respectively. 

As for the ROSAT PSPC, we selected a circular re- 
gion of 3'-radius around the source, while the data 
from an annulus of inner and outer radii 3' and 6' 
were used as the background. Since a super-soft source 
(SSS) IE 0056.8-7154 was located at 4' northwest of 
RX J0059. 2-7138 and the X-ray flux of this SSS was 
fairly large in the ROSAT band below ~ 0.5 keV, we ex- 
cluded a circular region of a 3'-radius around this SSS 
from the analyses of the ROSAT PSPC. 

3.1. Timing Analysis 

In total, the numbers of counts from the source regions 
are 21847 in the 0.7-10.0 keV band (ASCA GIS 2) and 
23250 in the 0.07-2.4 keV band (ROSAT PSPC). After 
the barycentric correction of arrival time for each event, 
we carried out timing analyses. 

The pulse period was already determined to be 
2.7632±0.0002 s with the ROSAT data alone (Hughes 
1994). We used the ASCA data to obtain a more precise 
period, because the duration of the ASCA observation 
is nearly five-times longer than the ROSAT observation. 
We performed epoch-folding on all of the ASCA GIS 2 
data, and determined a trial pulse period to be 2.76322 s. 
We then divided the ASCA observation into three seg- 
ments, each having a duration of ~ 33 ks, and folded 
each of the light curves with a trial period of 2.76322 s. 

To the folded pulse profiles of the first and the last 
segments, we made a cross correlation, and determined 
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the average apparent barycentric pulse period during the 
observation to be 2.763221±0.000004 s. This value is con- 
sistent with, and more precise than, the previous result 
of 2.7632±0.0002 s (Hughes 1994). 

The cross correlation between the middle and the 
first/last segments shows the upper limit of the phase 
difference to be ~ 0.04. Based on the assumption of a 
constant period change during the observation, this up- 
per limit of the phase difference is converted to the upper 
limit of the period derivative of \p\ < 6 x 10~ 10 s s _1 . 

Figure 2 shows the pulse profiles in the six energy 
bands (0.07-0.4 keV, 0.4-1.0 keV, and 1.0-2.4 keV in 
ROSAT PSPC, and 0.7-2.0 keV, 2.0-4.0 keV, and 4.0- 
10.0 keV in ASCA GIS) folded with a pulse period of 
2.763221 s. The pulse fractions, defined as (I max — 
^min)/2/avorago, in these energy bands are 0.07±0.06, 
0.14±0.04, 0.26±0.04 (0.07-0.4 keV, 0.4-1.0 keV, and 
1.0-2.4 keV in ROSAT PSPC) and 0.24±0.05, 0.37±0.05, 
and 0.36±0.05 (0.7-2.0 keV, 2.0-4.0 keV, and 4.0- 
10.0 keV in ASCA GIS), respectively. The pulse shape 
is nearly sinusoidal and independent of the X-ray energy 
The apparent decrease of the pulse fraction with decreas- 
ing energy, which was already found in the ROSAT data 
(Hughes 1994), is also seen in the ASCA data below 
- 2 keV. 

On a longer timescale (50-100 s), Hughes (1994) re- 
ported a flickering variability of about 30% of the mean 
rate in the ROSAT data. Such variability was also seen 
in the ASCA data with no energy dependence, while no 
larger variability, like a burst, was detected. 



3. 2. Spectral Analysis 

Figure 3 shows the phase-averaged spectra after back- 
ground subtraction. Since the pulse profile is essentially 
energy independent while the pulse fraction increases 
with energy below ~ 2 keV, we can naturally assume that 
the spectrum of RX J0059.2— 7138 consists of two com- 
ponents: a hard component dominant above ~ 2 keV, 
which has an energy-independent pulsation, and a soft 
component dominant below ~ 2 keV with no pulsation. 

With this working hypothesis, we first determined 
the hard component using the ASCA spectrum above 
~ 2 keV. Since no significant emission line was detected, 
we fitted the spectrum with a single power-law model. 
However, because systematic negative residuals remained 
above ~ 7 keV, we introduced a high-energy exponential 
cut-off to the single power-law model (ECUT power-law) 
and fitted the spectrum again. We then obtained an ac- 
ceptable fit to an ECUT power-law model with a photon 
index of T ~ 0.4, a cut-off energy of E c ~ 6.6 keV and a 
folding energy of Ef ~ 8.8 keV. 

To investigate the soft component, we then simulta- 
neously fitted the ROSAT and ASCA spectra in the 
0.1-10.0 keV band with various two-component mod- 




Fig. 2.. Normalized intensities of the folded light curves. The 
background levels are ~ 1% and ~ 10% of the average, in the 
PSPC data and in the GIS 2 data, respectively. 



els, in which an ECUT power-law model was always 
adopted to the hard component. For the additional soft 
component, we adopted a blackbody, disk blackbody, 
bremsstrahlung, thin-thermal plasma (Mewe et al. 1985), 
and broken power-law model. In these two-component 
model fits, a common absorbing gas column with solar 
abundance was adopted for both the hard and soft com- 
ponents. All models were statistically rejected at the 90% 
confidence level with wavy residuals below ~ 1 keV. In 
particular, an edge-like structure was seen at ~ 0.5 keV, 
which may correspond to the absorption edge of neutral 
oxygen atoms. Allowing the oxygen abundance in the 
absorption column to be free, we obtained statistically 
acceptable fits only when the soft component was either 
a thin-thermal plasma or a broken power-law. The best- 
fit parameters for these models are separately listed in 
tables 1 and 2, while the best-fit model is shown in figure 
3. 

We separately made phase-resolved spectra from 
phases 0-0.5 ('on-pulse') and 0.5-1 ('off-pulse'). We fit 
these spectra with the two accepted models, and found 
that only the normalization of the hard component varies 
with the pulse phase. The best-fit parameters of these 
models are also listed in tables 1 and 2, while the best- 
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Fig. 3.. Phase-averaged spectra. (Upper panel): Data are shown 
by crosses. The best-fit model (a thin-thermal plasma plus an 
ECUT power-law) of ROSAT PSPC and ASCA GIS are plotted 
with lines. (Lower panel): Residuals from the best-fit model. 
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Fig. 4.. Phase-resolved spectra. (Upper panel): Data are shown 
by crosses. The best-fit model (a thin-thermal plasma plus an 
ECUT power-law) of 'on-pulse' and 'off-pulse' are plotted with 
black and gray lines, respectively. (Lower panel): Residuals 
from the best-fit model. 

fit model is shown in figure 4. We also extracted the 
"pulsed spectrum" by subtracting the off-pulse spectrum 
from that of on-pulse. We found that the pulsed spec- 
trum is fitted with an ECUT power-law with absorption, 
of which the best-fit parameters are consistent with those 
of the hard component of the phase-averaged spectrum. 
These facts indicate that the pulsating X-rays contribute 
only to the hard component, consistent with the initial 
working hypothesis. 

4. Discussion 

Comparison to the Previous Results 
Hughes (1994) analyzed only the ROSAT data, and 
showed that the spectrum is fitted by a model having 



an extremely soft spectrum: a blackbody (kT ~ 35 eV) 
plus a steep power-law (r ~ 2.4). However, we have 
shown that the spectrum in the ASCA band (the hard 
component) is as hard as those of usual XBPs, and the 
ROSAT spectrum consists of mostly the "soft excess" 
below ~ 2 keV. 

For a consistency check, we first fitted our ROSAT 
spectrum with Hughes' model, and obtained best-fit pa- 
rameters consistent with Hughes (1994). However, as can 
be also seen in figure 2 in Hughes (1994), systematic pos- 
itive residuals were found above ^1.5 keV, which implies 
the existence of a hard component. We then fitted the 
ASCA and ROSAT spectra simultaneously with a three- 
component model, consisting of the hard ECUT power- 
law to the Hughes' two-component model. In this model, 
the best-fit blackbody temperature became much lower 
than the original result. After corrections for absorption 
and detector efficiency, the bolometric luminosity of the 
soft black body has an unrealistic value of ~ 10 45 erg s _1 . 

On the other hand, as we found, the combined ASCA 
and ROSAT spectra can be fitted with two different 
components, with a more reasonable luminosity of ~ 
10 38 erg s _1 for both components. Also, the energy- 
dependent pulse fraction and energy-independent pulse 
profile arc naturally explained by our two-component 
model. We thus infer that the present model is more 
probable, hence the spectrum of RX J0059.2— 7138 is not 
unusual for an XBP, at least in the hard energy band. 

4-2. Oxygen Over- Abundance 

We found an edge-like structure at 0.5 keV, which 
would be a hint of over-abundance of oxygen in the 
absorbing matter. One may argue, however, that the 
ROSAT spectrum fitted by Hughes (1994) does not show 
any structure at 0.5 keV. We should note that the black- 
body and power-law components in the Hughes's model 
cross at around 0.5 keV, which may produce an artifi- 
cial dip, and hence would compensate for the edge-like 
structure. 

To check whether the structure is real or not, we fitted 
a narrow-band spectrum (0.2-0.8 keV) with a power-law 
absorbed by cool matter of solar abundance. Because 
the structure remained near 0.5 keV, we claim that the 
edge-like structure is real. 

Since the oxygen abundance of the SMC interstellar 
matter is less than the solar value (Russell, Dopita 1992), 
we can infer that there exists oxygen-rich matter around 
the binary system of RX J0059.2— 7138. The companion 
star, which should have ejected the oxygen-rich matter, 
may be in an evolved stage where the stellar surface be- 
comes abundant only in oxygen, although such stars (e.g. 
WO stars) are very rare. 

On the other hand, metals such as carbon, nitrogen, 
and oxygen could be brought to the surface of a massive 
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Table 1. Best- fit parameters for a thin-thermal plasma plus an ECUT power-law model. 





kT 
(keV) 


z* 

(10 -2 solar) 


r 


E c 
(keV) 


(keV) 


Nl 
(10 20 cm" 2 ) 


z o 
(solar) 


4s 4h 

(10 3S erg s _1 ) 


X 2 /d.o.f. 


Av. 
On 
Off 


n 07+0.02 

u -°' -0.04 

o.35tg:8S 
o.38ig : g 


1 r.+0.8 
x -°-0.6 

1 6+ 1 ' 4 

1 fi +1 ' 4 


0.43±S;S 
0.4818$ 
0.47t°;}« 


6-5±i;i 

5-»±i:g 


9-3^ 
8.3™ 


4 3+0-3 
4 n+ - 3 

+ - u -0.3 


7 q+1.9 

6-4^1 

8.7il;| 


0.8 1.8 
0.8 2.2 
0.8 1.2 


63.9/62 
48.6/62 
59.4/59 



Note — The errors are at the 90% confidence level. 

*: Metal abundance of the thin-thermal plasma in unit of 1% of the solar abundance. 

f: Absorption column density with the solar abundance (except oxygen) in unit of 10 20 H cm -2 . 

\: Oxygen abundance in the absorbing matter in unit of the solar abundance. 

§: Absorption corrected luminosities of the soft and hard components in the 0.1-10.0 keV band at a distance of 
60 kpc, in unit of 10 38 erg s _1 . 



Table 2. Best-fit parameters for a broken power-law plus an ECUT power- law model. 



v\ r;* E* hl r 2 e c e { < z* 4 S 4, xVd-o.f. 









(keV) 




(keV) 


(keV) 


(10 20 cm- 2 ) 


(solar) 


(10 38 


erg s 1 ) 




Av. 
On 
Off 


z - D -0.7 

9 c+0.5 
z -°-1.5 
2 fi +0.4 


c 1+0.9 
°- 1 -0.6 
c 0+3.3 
°-°-0.9 

5 2+ 1 ' 8 


o.9±8:i 
o.9i8:? 
o.9l8:? 


n 4i + 006 

u - 4i -0.07 

46+ - 07 

u - w -0.08 

46+ - 12 

u -^°-0.16 


6.418:8 
6.511:? 

5-9±l:S 


9.0l 4 4 ° 
8.4liV 
8.1™ 


5.01S-J 
4-9l 2 i 
5.01S-J 


4 3+ 2 - 5 

4 3+ 4 -9 
^•°-2.9 
4 r+2.1 
4 -°-1.9 


0.8 
0.9 
1.1 


1.8 
2.2 
1.3 


63.0/61 
48.7/61 
57.4/58 



Note — The errors are at 90% confidence level. 

*: r and V are photon indices below and above E^, respectively. 

t: Absorption column density with the solar abundance (except oxygen) in unit of 10 20 H cm -2 . 
p. Oxygen abundance of the absorbing matter in unit of the solar abundance. 

§: Absorption corrected luminosities of the soft and hard components in the 0.1-10.0 keV band at a distance of 
60 kpc, in unit of 10 38 erg s _1 . 



star by turbulent diffusion, due to the rapid rotation or 
tides in a binary system (Maeder 1987). We thus fitted 
the spectra again, based on the assumption that abun- 
dances of carbon and nitrogen in the absorbing matter 
are the same as that of oxygen. The abundance was then 
determined to be ~ 6 solar, which leads us to a more 
comfortable possibility that the surface of the compan- 
ion star is enriched in CNO, which may be caused by 
tidal effects in the close binary system. 

We found that a thin-thermal plasma with very low 
abundance of 0.02 solar can fit the soft X-ray spectrum. 
As we discuss in the next subsection, the soft compo- 
nent is likely to originate from a large region surround- 
ing the binary system. This creates a dilemma that the 
X-ray emitting circumstellar medium is extremely under- 
abundant, while the X-ray absorbing circumstellar gas is 
CNO over-abundant. Although we have no clear idea 
of how to solve this dilemma, we suggest that the thin- 
thermal plasma is not a physical model, but is a phe- 



nomenological model like a broken power-law for the soft 
component. 

4-3. Origin of the Hard and Soft Components 

Since the pulsation was found to be mostly due to the 
hard component, we can conclude that the hard compo- 
nent originates from a small region, probably near to the 
polar cap of the neutron star. Hence, the origin would be 
the same as those of other XBPs. However, the cut-off 
energy is lower than that of the usual XBPs (10-20 keV; 
e.g. Mihara 1995a). According to the correlation between 
the cut-off energy and magnetic field strength (Mihara 
1995a, 1995b), the magnetic field of RX J0059. 2-7138 is 
estimated to be 6 x 10 11 < B < 1 x 10 12 G, which is one 
of the weakest magnetic fields in XBPs. 

A puzzle is the origin of the soft component. Since 
it shows no pulsation, it is likely that the soft com- 
ponent originates from a relatively large region, cover- 
ing a significant fraction of the full binary system. In 
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fact, when we adopt the thin-thermal plasma model as 
the soft component, the emission measure is very large: 
n 2 V ~ 10 61 cm~ 3 (table 1), where n and V are number 
density and volume of the plasma. If we assume that 
the plasma is distributed spherically, its radius should be 
larger than 0.07 AU under the optically thin condition, 
which does not conflict with the thin thermal scenario 
for the soft component. 

At present, only a few XBPs are known to have a soft 
excess. Her X-l has a soft component which is described 
by a blackbody (kT ~ 0.09 kcV; Dal Fiumc ct al. 1998), 
and is pulsating. McCray et al. (1982) found the phase 
shift between pulse profiles of the soft and the hard com- 
ponents, and argues that the soft component is produced 
by reprocessing of the hard component in the inner accre- 
tion disk. Another XBP, LMC X-4 has a soft component 
of a blackbody (kT ~ 0.03 keV) plus a bremsstrahlung 
(kT ~ 0.4 keV) (Woo et al. 1996). The blackbody com- 
ponent is almost constant during the pulse phase, while 
the bremsstrahlung component is pulsating with some 
phase shift from the hard component. Woo et al. (1996) 
argued that the blackbody component is emitted from 
the accretion disk, and the bremsstrahlung component is 
from somewhere near the neutron star, collimated as a 
fan-beam. 

The soft component of RX J0059.2— 7138 is, however, 
different from these XBPs; it is not pulsating, it is not a 
blackbody spectrum and it has no emission line. Hence, 
the origin of the soft X-rays would be different from those 
proposed by the above-mentioned authors. As a matter 
of fact, RX J0059.2— 7138 is a unique source among the 
XBPs with soft excess; because it probably has a Be star 
companion (Southwell, Charles 1996), it is different from 
either low-mass companion systems (Her X-l) or from 
high-mass supcrgiant companion systems (LMC X-4). 

Broad-band spectroscopic studies on the spectra of 
XBPs with a Be star companion have until now been 
very limited, due to the transient nature of this class 
and the relatively narrow energy ranges of the previous 
and contemporaneous satellites. New-generation satel- 
lites having a broad-band sensitivity and a wide field of 
view, such as ASTRO-E, XMM, and Chandra Observa- 
tory, will greatly improve our understanding on the soft 
excess of XBPs. 
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